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Introduction. The process of formation of fatigue defects in 
metal alloys with different structural morphology is consid- 
ered. The work objective is to develop a computational tool for 
determining the moment of the defect nucleation under cyclic 
loading. 

Materials and Methods. A physical model is built, calculation 
expressions are presented. The physical model is based on the 
theory of dislocations. It is shown that a structure factor is 
particularly important in the process of fracture nucleus origi- 
nation under dynamic cyclic loading. Depending on the struc- 
ture and properties of the material, as well as on the nature of 
the loads, the critical fatigue defect develops in the form of 
cracks, pores or micro-crater wear. 

Research Results. A numerical experiment was performed to 
determine the moment of nucleation of the critical-size defect 
in iron-base alloys under the drop hypervelocity impacts. 
Comparative data of calculations and bench tests for droplet 
impingement erosion of steels and alloys with the structure of 
ferrite, austenite, sorbitol and martensite are presented. The 
efficiency of the nucleation stage during the incubation period 
of erosive wear of the materials studied was evaluated. 
Discussion and Conclusions. There are no strict instrumental 
methods for determining the duration of the nucleation stage; 
therefore, it is recommended to use the proposed analytical 
model. In addition, the work performed gave a significant 
application result, i.e. it showed that the focused design of the 
material structure can significantly increase the wear re- 


sistance. 


Keywords: iron-base alloys, alloy substructure, cyclic load- 
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JloHCKOM rocyyjapCcTBeHHbIM TeEXHHYeCKHH yHUBepcuTeT, PocToB-Ha-Jlony, Poccuiickaa Deyepamna 


Beedenue. PaccmoTpeH tpolecc POPMUpOBaHHA YCTaJIOCTHBIX 
jedeKTOB B MeTAJIINYeCKHX CiaBax C pa3/IMYHOHW CTpyKTyp- 
HOH MopdonormeH. Llenb padoTbl — co3jjJaHve pacueTHoro 
alnapata Jif OllpezesIeHHA MOMECHTA 3apOxK CHA yKa3aHHBbIx 
TeeKTOB B YCJIOBHAX IMKIIM4eCKOLO HarpyxKeHHuaA. 
Mamepuaibi u memoovi. IloctpoeHa du3u4eckad MOJelIb, 
IIpeJCTaBJICHbI pacueTHble BbIpaxkeHua. B ocHoBy ~u3su4ecKoOn 
MOeIM MOOKeHAa Teopua AHcoOKayHH. Iloka3aHo, 4To mpu 
MHAMWYeCKHX WMKIMYeCKHX HarpyKeHHAX OMpeesArolsee 
3HadeHHe B Ipollecce 3apoxeHHA O4ArOB pa3pylleHnaA UMeeT 
CTpyKTypHbIH (daxTop. B 3aBHCHMOCTH OT CTPyKTypbI U 
CBOHCTB MaTepualla, a TaKxKe OT XapakTepa Harpy30K, KpHTH- 
YeCKHM YCTAaIOCTHbIM edeKT pa3sBuBaeTcA B (POPMe TPeLIHHBI, 
MOpbl WIM MUKpOKpatTepa H3HOCAa. 

Pe3yibmamol ucciedoeaHua. BbinosIHeH 4HCIICHHbIM 9KCTIe- 
PUMeHT TO ONpeyeueHHIO MOMeHTa 3apoxKTeHHA WedeKta 
KPHTHYeCKOTO pa3Mepa B CliyIaBaxX Ha OCHOBe >%*esle3a Mpu 
BbICOKOCKOPOCTHBIX KalleJIbHBIX COyqapeHuax. ITpeacTaBieHsl 
CPaBHUTeJIbHble JAaHHble pacueTOB MW CTCHJOBbIX UMCIIbITaHnM 
10 KallleyqapHou 9po3HU CTasleH MU CIIaBOB CO CTpyKTypou 
(bepputa, aycTeHuta, copOutTa u MapTeHcuTa. OleHeH BKIaI 
CTaJIMU 3apoOx*K CHA B MHKyOalMOHHbIN NepvO 2pO3HOHHOTO 
W3HOCAa MCCIEOBAHHbIX MATepHasIOB. 

O6cyocdeHue u 3aKlOUeHUe. CTporue WHCTPyMeHTAaJIbHbIe 
MCTOJIbI JIA OMpeeICHHUA MpOOJDKUTeIbHOCTH CTaqWu 3a- 
POX CHUA OTCYTCTBYIOT, IOITOMY PCKOMCHOBAHO UCIIOIIb30- 
BaTb I[peJIOX%KeCHHYIO paC4eTHO-aHasIMTHYeCKy!O MOJICIIb. 
Kpome Toro, BbIMOJHeHHad padota Jayla BaxKHBIN NpukTaqHoU 
pe3ylIbTaT — IIpOeCMOHCTpupoBala, YTO WeseHallpaBsIeHHoe 
KOHCTpyHpOBaHHe CTPyKTypbI MaTepHasia MOXKeT CylIeCTBeH- 


HO YBeCJIM4YHTb U3HOCOCTONKOCTB. 


KuroueBble CJ10Ba: CIiIaBbI Ha OCHOBe 2KesIe3a, CyOCTpyKTypa 
CIlaBOB, IMKIMYeCKHe Harpy3KH, 3apoxKTeHHe O4AaTOB pa3- 
pylieHus, (u3n4eckad MOj[esIb, UYHCICHHBIM 3KCIICpHMeHT, 
KaluleyapHad 9po3H4. 
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Introduction. Modeling specific unpredictable phenomena is a pressing challenge in terms of preventing their 
undesirable consequences. For example, in aeronautics, heat power engineering, gas turbine production and agriculture, 
the problem of droplet impingement erosion remains unsolved. The issue is discussed on the world level [1—4]. The first 
success in modeling this phenomenon concerned single0020 droplet collisions [5—10]. Some papers [11—13] attempted 
to quantify the ability of materials and coatings to resist the action of dynamic cyclic loads. In particular, it was possible 
to create an analytical model tested under the conditions of liquid droplet shock [14—18]. Many authors [19-21] do not 
differentiate the behaviour of a material or coating at the stages of nucleation and development of fracture. Nucleation is 
problematically identified in the experiment. In the model mentioned above, it was estimated using empirical coeffi- 
cients that did not have a universal theoretical foundation, and this is a clear drawback of the solution. This work objec- 
tive is a deeper theoretical study of the stage of the defect nucleation in the coating causing its wear and destruction, as 
well as the creation on this basis of a calculating apparatus for determining the moment of nucleation of cracks or pores 
in various materials and coatings under cyclic loading. 

Materials and Methods. Difficulties in evaluating the stage of fracture and wear nucleation under the repeat- 
ed exposure (for example, at high-speed liquid-droplet collisions that provoke droplet impingement erosion of metal 
products) are associated not only with instrumental determination of nuclei of cracks or pores. At this stage, the pres- 
ence of structural elements in the material does not allow us to present the process as a continuum and to use linear frac- 
ture mechanics (for example, the Paris-Erdogan fatigue theory) for calculations. Therefore, in the fracture mechanics, 
the initial stage of the fatigue defect nucleation is singled out in a special stage. It is called “microstructural fracture 
mechanics” and proceeds until the defect reaches the size of /, = (4...10)-d, where d is the size of a structural element of 
the material [22]. An analytical estimation of the stage of fatigue defect nucleation requires, first of all, a metal physical 
approach. The task is to find a calculated expression for the number of cyclic loads N3 necessary for the formation of a 
critical size fatigue defect in the material. As an empirical foundation for the computational model, the experimental 
data of bench tests for droplet impingement erosion, where the number of drop collisions plays the role of N3, can be 
used. 

In the general case, the equation for the number of dynamic cyclic loads N3 should include three components: 
mechanical, kinetic, and structural. 

The metal physical concept of a mechanical component is related to the number (density) of mobile disloca- 
tions p,, arising from a single shock loading (collision). 

. In the process of multiple collisions, the number of mobile dislocations increases, they move in the metal ma- 
trix along slip planes under the action of tangential stresses o, until they form flat clusters with critical density p,, at the 
nearest impassable barriers. Exceeding p;,. causes spontaneous breaking of interatomic bonds in the metal and the origi- 
nation of a crack nucleus. The p;,, values for various materials and coatings are known [23]. A mechanical component of 
the number of collisions is expressed in the form of the dependence N3 = f (pyx-/P mn). 

A kinetic component of the number of collisions V3 takes into account the dependence of the dislocation ener- 
gy on its velocity V,. From the theory of dislocations, it is known that with an increase in speed of dislocation, its ener- 
gy increases according to the expression of Einstein for bodies moving at speeds close to the light speed. The only limit 
for the dislocation velocity is the sound speed in a C) crystal, at which the dislocation energy becomes infinite. There- 
fore, considering the specifics of the drop impact, a kinetic component of the number of collisions is expressed as 
N3=f (Val Co). 

A structural component of the number of collisions NV; has two aspects. The first is that the dislocation motion 
under impact is limited by the size of the structural element of the D matrix, within which a free path of dislocations is 
possible. In the general case, a grain size is considered as such an element, and its thickness — for thin coatings. The 
second aspect considers that within the grain or coating, there can be obstacles to the dislocation motion: second-phase 


particles, small-angle boundaries, immobile dislocations fixed by atmospheres, stacking faults, lattice resistance (Peierls 
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stress). As noted above, these both aspects are reflected in the structural component of the number of collisions, and it is 
a function of two variables N3 = f(D, AG), where AG is the Gibbs free energy to activate the process of obstacle negoti- 
ation within the structural element D. The negotiation mechanism depends on the nature of the obstacle, which deter- 
mines the value of AG for each specific case of the material structure (their classification is given in [23]). 

Thus, all the components of the desired number of collisions N3 affect the process of fatigue defect nucleation 
at the same time — with each collision. This fact determines the commutative nature of the interaction of the mechan1- 


cal, kinetic, and structural components. Then, in the most general form, N3 1s determined from the expression 


2 

V Dp. a 

N, =F. {2 ae (1) 
Pn Cy) 


where & is the Boltzmann constant; Tis thermodynamic temperature, K; /p is the path covered by mobile dislocations in 
one loading cycle (collision). 

The expression (1) is obtained on the basis of the theory of dislocations and represents a theoretical concept. 
To use (1) in assessing the wear resistance of materials and coatings under the conditions of cyclic liquid-droplet colli- 
sions, it 1s necessary to specify its constituent quantities (such as p,,, Va, /), AG) through measured impact parameters 
(e.g., impact velocity Vo and droplet size Ry). The authors have already completed the corresponding refining of the 
model (1) and are preparing material for publication. 

Research Results. The computational model has been tested for iron-based alloys of various compositions and 
structures. This selection is due, first of all, to the difference in the behaviour of dislocations in these alloys and makes 
it possible to identify the capabilities of the computational model. So, in ferrite and austenite, the factor hindering dislo- 
cations is the lattice resistance (Peierls stress). In a sorbitol structure, the glide of dislocations is determined by the mor- 
phology and distribution of obstacles — by carbides Fe and Cr. In martensite, the dislocation motion is hindered as giv- 
en. Thus, the subject of research in this part of the work is actually a structure factor. 

Since all the materials studied are alloys based on iron, their basic stress-strain properties differ insignificantly. 
The following values were used in the calculations: elastic modulus EF = 186 ... 218 GPa; Poisson ratio v= 0.20 ... 0.31; 
shear modulus 1 = 64 ... 80 GPa; Burgers vector b = 2.5-10'° m. The data of numerical experiments on the implementa- 
tion of the presented model and bench tests of Fe-based alloys with different structural morphology, including hetero- 
geneous ones, are given in Table | and in Fig. 1. 

Table 1 
Experimental mo and calculated N3 values of the number of collisions* for Fe-based alloys 


Collision parameters 
Structure parameters Vy = 250 m/s; Vy) = 340 m/s; 


Material Ry = 0.55 mm Ro= 0.32 mm 


ie a ja, |ptm'y| rims | eee} Pap _Eateon sr 
i me te 


ARMKO F 
~ | BCC (2.9) 0.5 100 3.0 2197 1187 
(ferrite) 


O8X18HI10T 
FCC (3.06) pade || ames | Seen) 255s || basa | S55 
(austenite) 
20X13 0.35 
| Bcc (2.9) | 2.0 | 100 4655 | 0.280 | 16630 | 3742 | 0.254 | 14768 
(sorbitol) jpocias)| 20 | 100 | x, | aoss | 0280 | rosso | ara | o2ss | 14768, 
20X13 Teas. 
| eee 0.02 | 0.01 0.01 9344 | 0.550 | 17861 | 5705 | 0.470 | 12140 
(martensite) (2.95) 


* Number values of droplet collisions falling at one point are presented (determined in the software of the test bench). 
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c) 
Fig. 1. Dependences of the number of droplet impacts N3 necessary for crack nucleation as a function of impact rate Vo at Ro = 0.55 
mm for Fe-based alloys with different structures: a) ARMKO Fe (ferrite); b) steel sorbitol 20X13; c) steel martensite 20X13 


Values of some parameters given in Table | need explanation. The distance between obstacles / on the path of 
mobile dislocations in ferrite and austenite is the distance between dislocations, 1.e. / = 1/ Jp , - Given high plasticity of 


these solid solutions, the density of their dislocations is relatively low pz = 10'' m°, which gives /~ 3 um. In hardened 
steel, the dislocation density is close to critical p,,= 1016 m~, which not only gives the structural parameter value 
/= 0.01 um in Table 1, but also changes the design procedure for the number of collisions in the source model (1). 

For hardened metal materials (for example, for strongly deformed or hardened alloys including 20X13 with a 
martensitic structure), in the expression (1), the first three factors become a unit. In such materials, the critical disloca- 
tion density p,,. has already been reached; therefore, the concept of mobile dislocations p,, loses its meaning. In terms of 


numbers, this is expressed as follows: px/Pm=1; Vi=0; D=1, = 1/ - Ip, , 


As for the exponent in the expression (1), it retains its meaning with the values of the parameters a, and / giv- 
en in Table 1. A crack nucleates in such a structure, an external cyclic action should overcome the existing stresses in a 
saturated dislocation medium and form a shear defect. Given the above, the designed expression (1) for hardened steel 


takes the form [23-26]: 
oY Ae | gee 
N, = ag J al 5) 0) 


Here, AF is the activation energy of the process of overcoming obstacles without applying external voltage, the value 
AF actually determines the strength of the obstacles from the point of view of their dislocation; o, is the internal stress 
existing in the material and allowing dislocations to pass through the obstacle at minimum AG values. The yield stress 
of a solid at T= OK 1s taken as o,. The values o, and AF are the material properties and, in the general case, are ex- 
pressed through its basic physical characteristics u and b: 
Go, =p: b/EAF=aq,- ud. (3) 
In the expressions (3), the coefficient a, classifies obstacles according to their strength [23]. 
In the expression (2), the shear stress in the slip plane under the action of which the mobile dislocations move, 
is considered as the stress o,. Stress o, is the projection of the external stress vector onto the slip plane and is related to 
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the lattice slip system through the Taylor factor M, according to the Schmidt - Boas law for polycrystal [27]. It can be 
expressed in terms of impact parameters Vp and Ro [28]. 

The experimental values shown in Table | characterize the incubation period of erosive wear of the material 
mo. This value includes both the stage of the defect (cracks, pores, microcraters) nucleation, and the stage of its devel- 
opment before the sample loss-of-mass start. It does not seem possible to instrumentally identify the number of colli- 
sions occurring only at the nucleation stage in the experimental values of mg; therefore, the values of N3 are exclusively 
theoretical. The contribution of the nucleation stage to the total value is characterized by the coefficient a») = N3/ mo. 

In addition to the numerical values presented in Table 1, the model considers the number of collisions N3 nec- 
essary for the fracture nucleus origination (cracks, pores, microcraters) and the impact rate Vp at the specified droplet 
size Ro. In addition, using this solution, we can determine the dependence of V3 on Vo. Some options of such depend- 
ences are presented in Fig. 1. 

The results show that the designed data N3 agree with the experiment (the bench test data mo). They also 
demonstrate compliance with the basic canons of fracture mechanics. In particular, the coefficient a» characterizes the 
nucleation energy with respect to the fracture energy. In plastic metallic materials (ARMCO, austenite, sorbitol), the 
nucleation energy of fatigue cracks is significantly less than the energy of their growth. In hardened alloys (martensite), 
the nucleation energy is almost always higher than the growth energy. As a rule, this ratio is much larger in favour of 
the nucleation stage. Table 1 data show that in steel 20X13 with a martensitic structure, the stage of nucleation of fa- 
tigue defect N3 is about half of the entire incubation period of the wear formation mg: Q») = 0.47 and 0.55 for collision 
velocities Vo = 340 m/s and 250 m/s, respectively. Whereas for the same steel with a sorbitol structure, the coefficient 
Qo 1S half as much. The values of a obtained by us correspond to the data of other authors for similar conditions of cy- 
clic loading (for example, for steel 30XI‘CH2A in [29]. 

The computational dependences in Fig. 1 also correspond to the experimental data. They are located asymptot- 
ically with respect to the values Vy = 100 ... 150 m/s (the numerical experiment did not study this region in detail). The 
computational dependences mentioned above correspond to the results of bench tests, which showed that at 
V)< 135 m/s, wear is generally not observed in the iron-based alloys [30-34]. That is, the asymptotic property of the 
graphs in Fig. 1 confirms this empirical fact: at the indicated low collision velocities, fatigue defects of a critical size do 
not nucleate. 

Discussion and Conclusions. There are no instrumental methods for accurate determination of the stage dura- 
tion of fracture nucleus origination of materials, therefore, if necessary, the proposed analytical model is used. It is ap- 
plicable to materials of various morphologies, in which the mechanisms of defect nucleation (cracks, pores, microcra- 
ters) have a dislocation nature. Therefore, under normal conditions, the model is not applicable to ceramic materials 
with a high proportion of covalent bonds. Another use restriction for the model is, perhaps, too high sensitivity of the 
exponential factor in the expressions (1) and (2). A small error in determining the numerical values of a, o, or / in the- 
se formulas prevents from obtaining the results adequate to the experimental data. 

The performed work yielded an important applied result, it showed that the focused design of the material 
structure can significantly increase wear resistance. 
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